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I. INTRODUCTION 



In this paper we concentrate on the few-body interactions of 77 mesons. These could 
complement our knowledge on the 7/-nucleon interaction and tell us possible evidence of 
//-nuclear quasibound states. Such quasibound states were predicted by Haider and Liu |T| 
and detailed calculations performed by Li et al. 0, when it turned out that the 77-nucleon 
interaction was attractive. To be observable these states should be narrow enough, and this 
is not likely to happen for the lowest 77 states in large nuclei. On the other hand it was 
suggested by Wilkin P that the rapid slope seen in the pd — > ?7 3 He amplitude of Ref. || 
just above the threshold may signal that a quasibound state is generated already for small 
nuclei(A = 3). In contrast, a recent study of the dd — > ?7 4 He reaction shows no structure 
due to any final state r/ 4 He interaction ||. All this could indicate a large ?7 3 He scattering 
length and a small one for r/ 4 He. However, quite an opposite interpretation is put forward 
in this paper. We calculate the r/ 3 He and 7/ 4 He scattering lengths and find that the former is 
smaller than the latter, and that they also differ in the sign of the real part. This suggests 
that the 77-nucleus attraction is not strong enough to give any binding effect in the ?7 3 He 
system, but it is likely to give one in the 77 He system. 

In the standard theory of final state interactions the energy dependence of reactions 
is assumed to be determined by the scattering amplitude between the final state particles 
||. In this paper we show that the shape of the low energy 77 production cross section is 
also significantly influenced by an interference of the free and the scattered waves in the 
final 77-helium states, because the corresponding scattering lengths are not very large. This 
interference is such that the decrease with energy becomes steeper for both 3 He and 4 He than 
that calculated from the final state scattering amplitude alone. However, in the scattering 
amplitude itself the real and imaginary parts of the scattering amplitude, due to the above 
mentioned difference in the sign of the real parts, could be expected to conspire so that the 
slope in the 4 He case would be somewhat smaller than for 3 He. Numerical results do not 
support this for r/N scattering lengths considered realistic. 
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Before this physical interpretation of final state interactions is discussed in Section |HI| , a 
formalism is developed in Section [II] to calculate the ?7-helium scattering lengths. By some 
formal manipulations the multiple scattering series is summed. The procedure used is shown 
to converge quickly in the case of r/He optical potentials, which may be solved exactly using 
the Schrodinger equation. Then, necessary corrections to the optical potential limit may be 
easily implemented by modifying the equivalent multiple scattering series. 

The conclusions are not fully quantified since the 7/-nucleon input is not determined 
uniquely. Also the i] production mechanism is not under full control. Here a method for 
calculating only the final state interaction is given. However, this method is presented in 
sufficient detail that a more complete comparison with the data and the determination of 
the input uncertainties can be made when measurements proposed at proton storage rings, 
such as the one at Celsius, are performed. Comparing specifically with the r] 3 He data we 
also perform an extensive variation of the r]N scattering length in search for a constraint on 
it, complementing elementary photo- or electroproduction. 



II. SCATTERING LENGTHS 

The ^-helium scattering lengths are calculated in this section. At first we consider the 
simplest optical model expressed in terms of the 77-nucleon scattering length. This problem 
may be solved numerically, but in order to improve the method an equivalent alternative for 
the optical model is provided, which consists of a partial summation of the multiple scattering 
series generated by the optical potential. The sum is expressed in terms of multiple integrals. 

Next, some necessary improvements to the optical model are introduced into the partial 
sum. These are essentially twofold: (i) removing multiple collisions on the same nucleon 
and, (ii) introducing an off-shell 7/-nucleon scattering matrix. Other effects, e.g. the Pauli 
principle, are less significant. These improvements introduce massive changes to the 77- helium 
scattering lengths defined here as the zero energy limit of the effective range expansion 

pcot5 = - + ^r p 2 . (1) 
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In particular, for the r] 3 B.e system, the simplest optical model yields a large negative real 
part of the scattering length indicating the existence of a quasibound state. However, when 
the corrections are included, a sizable positive scattering length emerges. On the other hand, 
in the r/ 4 He system we do find indications for a quasibound state close to threshold. 

A. Multiple scattering expansion for the inverse scattering length 

In Ref. a multiple scattering scheme was proposed to calculate the energy shifts and 
widths in the atomic states of antiprotons interacting with a light nucleus. In this paper 
we apply the same method to r/He scattering at threshold. First the procedure is presented 
in some detail, since it is important to understand to what extent the basic form of the 
multiple scattering scheme is, in this case, numerically equivalent to the standard optical 
model approximation. 

The scattering matrix T(i]A) for an 77 meson interacting with a nucleus of A nucleons 
may be expressed as a series in the following way by first considering the scattering from 
two non-overlapping fixed centres. In that limit it can be shown that the scattering matrix 
at zero energy has the exact form || 

T _ T, + T 2 

1 1 (T 1 +T 2 )D(T 1 +T 2 ) ' 
2 (Ti+Tb) 

where D — l/l is the propagator of the scattered particle - with I being the distance apart 
of the two scattering centres - and the Tj are the scattering matrices from the separate 
centres. This expression has the following feature that is important in few body systems. By 
expanding the denominator, a multiple scattering series emerges which - through the factor 
of 1/2 in the denominator - automatically takes into account the exclusion of successive 
scatterings from the same centre. Being guided by this and denoting T = J2Ti, one might 
then naively expect the analogous expression for the scattering from A non-fixed centres to 
have a form 

T(r)A) =<T> [l + P]-\ (3) 
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where 



(A-1)<TDT> 

A<T> {) 

with 

D = — m/(27r|r — r'|) (5) 

the zero-energy propagator for a free rj. Here the factor (A — 1)/A is the generalisation of 
the exclusion factor 1/2 in Eq. (0) and m is the reduced mass of the r\A system. 

However, as it now stands Eq. (|^) cannot be correct beyond 0(T 2 ), since it does not 
give the required form 

A-l\ „, /A-V 2 



T(r]A) =< T > + ( — — J < TL>T > + J— — J < TDTDT > 

+ (A^) < TDTDTDT > +.., (6) 

when expanded in powers of T. One way of ensuring that this correct expansion results is 
to modify Eq. to 

T( V A)=<T> [1 + P + Q + R + ..]- 1 , (7) 

where the quantities P,Q,R,.. are of order T, T 2 , T 3 , respectively, and are chosen in turn 
to guarantee Eq. @. For example, on expanding the denominator of Eq. (0) the term of 
0{T 2 ) is -Q + P 2 , which in Eq. (|) should give (^f < TDTDT >. Since P has already 
been fixed by the second term in Eq. we get Q = P 2 — < TDTDT >. This 

is a unique procedure and, neglecting for the moment the above exclusion of consecutive 
scatterings on the same nucleon, it leads to the expressions 

< TDT > < TDT > 2 < TDTDT > 

1 ^£ 



<T> ' ^ <T> 2 <T> 

< TDT > 3 < TDTDT >< TDT > < TDTDTDT > 
R ~~ <T>3 +2 ^TT 2 Tt^ • (8) 
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for the first three terms of the series in the denominator. Also in Q,R... it is immediately 
seen that the integers multiplying each term cancel each other e.g. in Q we see +1 and -1, in 
R we see -l,+2,-l etc. In the fixed centre limit all of the integrals reduce to the same value, 
so that Q,R... are then all zero as found in Ref. ||. The important point is that numerically 
this cancellation continues to a great extent even away from the fixed centre limit as seen 
below. Therefore, one could hope that the introduction of the ratios < TDT > / < T > 
etc. of double and triple scatterings and other " disconnected" terms at various places would 
speed up the convergence of the multiple scattering series. 

Here now T denotes a scattering matrix of the rj from A nucleons in the impulse approx- 
imation. At "zero" energy mT(f]A)/2ir reduces to minus the i]A scattering length a(rjA), 
and 

T = —tAp{r) (9) 
A 4 

with t being the ^-nucleon scattering matrix at the appropriate energy, Ap(r) the nuclear 
density and p the reduced mass for the rjN system. The expectation values appearing in 
Eqs. (H) can be expressed in terms of the propagator and nuclear density as 



<T>= —tA [ dr p(r) = —tA (10) 
p J p 



< 



TDT >= -| [jtJ^J J drdr'p(r) F ^p(r') (11) 



< TDTDT >= 2 ^(™t) 3 A3[ [ [drdrW' ^P^l , etc. (12) 
m \p J J J J | r — r 1 1 r — r" | 

Using the Gaussian density profile 

p(r) = l/(v^i?o) 3 exp[-(r/i? ) 2 ] (13) 
one obtains now the expansion coefficients 



P -t( Am \ Q -t 2 ( 



Am 
RrmsH J 







)• 


f--ll 




.7T 



R = t 



( Am V 



^ 2(-) + 0.7796^1-) 



(14) 



where -Rrms — y3/2i2 is the RMS matter radius of the A nucleons. The number 0.7796 
in R is the result of a double summation and is expected to have an accuracy of ±0.0001. 
When, the terms of the series in the denominator are clustered into increasing powers of t as 
indicated in Eqs. d^H), it is found that there exists a considerable amount of cancellation, 
e.g. P/d = 0.977, Q/d 2 =-0.045 and i?/d 3 =0.0076 where d = t (Am/ R RMS p). Therefore, if 
t/RmAS is reasonably small - as it is in the present case of rjN scattering - the series appears 
to converge rapidly. 

A check on the convergence is given in Table H, where a comparison is made between 
the above series expansion for T(r)A) and its value calculated directly from the equivalent 
optical model potential 

2-7T 

V(opt) = — t ■ A p(r) . (15) 

In this comparison the factors (A — 1)/A in Eq. (|5|) must be neglected. For completeness, 
the impulse approximation (IA) result 

Tfi 

a(r}A, I A) = —A — t (16) 
A 4 

is also quoted in Table | [i.e. T(rjA) with P — Q = R = 0]. As a first approximation t is 
taken to be minus the rjN scattering length, i.e. t(E = 0) = —a(rjN). The actual numbers 
used are a representative sample from the following sources. Several groups have performed 
coupled channel analyses of ^-nucleon and 7r-nucleon scattering [l|-^|,p-|TT| . These differ in 
the input data and also in some details of the extraction of the a(rjN). The rjN scattering 



lengths obtained are: 0.27+i0.22 and 0.28+i0.19 and (0.50 ± 0.20) + z(0.33 ± 0.06) fm 



]. The recent electroproduction data yield 0.476+i0.279 fm [[]], while photoproduction 



experiments suggest the possibilities: 0.430+i0.394, 0.579+i0.399, 0.291+i0.360 fm |TT . 

The multiple scattering effect is dramatic as compared with the impulse approximation 
alone, changing the attractive real part of the rjN amplitude into a repulsive real part for 



the t]A amplitude. However, the near equality between a(r]A) and a(opt) gives confidence 
that, indeed, the series in Eq. (|7|) is rapidly convergent and the use of only the terms P, Q 
and R gives a sufficient accuracy. 

Having shown in some detail that, indeed, the standard optical model approach can be 
replaced by the multiple scattering series of Eq. ([?]), it now seems justified to modify the 
latter to include effects not so easily incorporated directly into the optical model. 



B. Corrections to the optical model 

Several improvements can be made to the optical model approach and these can be 
implemented into the partial sum of Eq. (|7|). 

i) Firstly, as shown in Eq. (§) the factors A 1,2,3 in the rescattering quantities P, Q and R 
calculated in Eqs. ([14]) should be replaced by (A — 1) 1,2 > 3 to prevent the rj from interacting 
successively with the same nucleon, i.e. a(r]A) — > aA-i{i]A). In Table |I| this effect is 
demonstrated for the rjN scattering lengths used in Table |, and it is seen to have a large 
effect in all cases. In particular, this correction makes the real parts of the scattering lengths 
small for 3 He making the existence of a quasibound state in the rj 3 He system indicated by 
the optical model questionable. The absence of such a state seems to be further confirmed 
by later corrections for 3 He. However, the real parts tend to become even more negative in 
the case of 4 He. 

ii) Another improvement to the above series is to use an rjN scattering amplitude that is 
more appropriate for scattering on a bound nucleon in a medium. This can be approximately 
taken into account by extrapolating a(rjN) off the energy shell through replacing the above 
scattering amplitude a(0) at zero energy by a(off) at a negative energy defined via the 
equation 

1 1 1 



a(0) a(off) a(0) 



iK„, (17) 

where K v = i^j2fi(E scp + E TCC ) with -E^rec being the A — > (A — 1) + 1 separation energy 
and the recoil energy of the rjN pair relative to the residual nucleus. For 3 He [ 4 He] these 
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quantities have the values E sep = 7 [21] MeV and E rec = 12 [12] MeV. The effects are shown 
in Table |TTT]. There it is seen that our best estimate aA-i{f]A, off) of the r]A scattering length 
is quite different from that predicted by the optical model. For example in r/ 3 He scattering, 
the negative Rea(opt) has turned positive and comparable to that of the original Impulse 
Approximation. This indicates that there is no binding in this system. However, for ^ 4 He 
the negative sign of Re a(opt) is maintained suggesting a quasibound state. At the end of this 
section these effects are interpreted in terms of poles in the scattering matrix. One should 
note in this context that, since in Eq. fll7|) a significant non-zero value is assigned for the rj 
momentum, the next term in the effective range expansion ([!]) could become important, if 
r is large. This could have now the effect of changing the energy variation present in Eq. 

(0)- 

iii) The major mechanism that generates the imaginary part of a(r]A) is the reaction 
rjAi — > N*(A — 1) — > TfAf, where N* is the nucleon resonance iV*(1535) with a strong 
coupling both to the rj and the pion. Therefore, for rj scattering on deuterium or 4 He - both 
isoscalars - the final nucleus Af cannot be an isoscalar. Because the spin is not involved in 
this s-wave scattering, then for example with the deuteron the final iViV state must be the 
3 Pi state and also the transition operator must be spatially antisymmetric. This opens up 
the interesting possibility that pionic inelastic channels are damped in these cases leading to 
a reduction of the in-medium value of Im a(rjN). However, as shown in the Appendix this 
turns out to be only a very small effect and so this correction is not included in the present 
calculations. 

If there exists a pole in the scattering matrix close to the threshold, the scattering lengths 
may become larger than the nuclear radius. To some extent this situation is met here, in 
particular in 4 He. In the case of a bound state Re a < 0, while a virtual state corresponds to 
Re a > 0. The connection is unique provided the effective range is small, which is assumed 
here. However, the validity of this assumption is not clear. Another complication arises 
because of the presence of decay channels described here by Im a(j]N). Even though there 
is no detailed many channel structure of the scattering matrix, let us, however, look for the 
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poles given by the condition (1 — ipa)— 0. With our best values aA-i(i]A, off) we have a pole 
in the upper complex momentum half plane, i.e. a quasibound state in the ^ 4 He case. This 
situation is, in fact, typical for all Re a(j]N) in the range 0.3 to 0.6 fm or even higher. On 
the other hand, with positive Re (Ia-\ (tjA, off ) one meets a virtual state in rj 3 He systems. 



III. FINAL STATE INTERACTIONS 

Since there are no beams of ^-mesons, the interactions of these mesons may be seen only 
via final state interactions or via the decay mechanisms of quasibound states. As seen in Fig. 
1, the pd —>■ i] 3 He production amplitude displays a rapid fall-off away from the threshold 
region, which led Wilkin to conjecture that an r] 3 He quasibound or resonance state exists 
nearby 0. This is reflected by the approximate proportionality of the cross section to the 
final state interaction factor || 

a{r)A) 



(18) 



1 — ipa(r]A) 

where a(r]A) is the 77- helium scattering length and p is the 77 momentum. It was found by 
Wilkin in the optical potential approach ||, recalculated here in Section O, that Im a(rfHe) 
is rather large, which gives the required slope and indicates a singularity. However, surpris- 
ingly the recent data on the reaction dd —>■ r/ 4 He indicate no such slope in the cross section 
close to the threshold ||. We now analyse these two measurements below. 

First, let us note that Eq. ([18]) provides a good description only, if |a(77A)| ^> -Rrms ~~ a 
condition not well satisfied here by the -Rrms for 3 He. A more general model needs the final 
s-state wave function for the 77- He system ip~(r). One particularly simple form of %p~(r) is 
that from a separable potential with the Yamaguchi form factors (1 +p 2 /j3 2 ) [18], which 
gives 

^.-^ _ sin(pr) | [exp(-ipr) -exp(-(3r)} 
pr ' r 

Here / = Fi — a(r]A)/(l — ipa(r)A)) is the on-shell ?7-helium scattering matrix, where the 

a(r)A) are taken to be the cia~i{i]A, off) from Section |T| and not the a(r]A) given by the 
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separable potential. Since the factor [exp(— ipr) — exp(— /3r)] determines the behaviour of 
the scattered wave inside the range of the interaction, it can be interpreted as producing an 
off-shell effect into the reaction. A plausible choice of (3 — 1/.Rrms is taken - but, as shown 
below, the shape of the cross section is rather insensitive to the actual value of (3. 

In the reaction process the 77's are produced with some amplitude H(r,pi) that depends 
both on the initial projectile momentum (pi) and on the spatial extent of the process. For 77 
energies in the range of 0-5 MeV the dependence on pi («1 GeV) is presumably small. So far 
there is no complete understanding of the actual production mechanism P,|T5l . However, for 
the present purposes it is sufficient to make only some rather general qualitative statements 
concerning this mechanism. Here we simply assume a proportionality of the production 
amplitude to the nuclear density used to derive Eq. ( |14j ) H{r) = exp (— (r/i? ) 2 ) with 
Rq = A^/2/3-Rrms an d A ~ 1 being a natural choice. In this way, the final state interaction 
factor becomes 



At first sight it appears that this model for incorporating final state interactions contains 
two adjustable parameters A and (3. However, in practice the f3 dependence is weak with 
even (3 = 00 being not unreasonable. As said above, we typically fix (3 at 1/Rrms leaving 
only the A dependence. Lacking an actual model for 77 production, in all cases the results 
are normalised to give the experimental value of the spin-averaged quantity 



at p v = 0.246fm . 

The original hope had been that, with A around unity, a good fit would be obtained to 
the shape of the experimental data. However, this was only so for potential IV. In that 
case, with A=0.88 and f3 = 1/Rrm.s there was a very shallow minimum in the \ 2 fit f° 
|/(expt)| 2 . It should be added that this fit did not include the lowest experimental point at 
Prj = 0.051fm _1 , since this is thought to be subject to large systematic errors due to beam 




(20) 



/(expt)| 2 



Pd da 

Pr, dfl 



(pd -> r? 3 He) = 0.63 ± 0.02 //b/sr 



(21) 
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width effects, including energy losses in the target || . The results are shown in Table p"V 
and Fig. 1. 

This table illustrates the following points: 

1) As seen from columns 3 and 4 the dependence on (3 is weak. Both (3 = 1/Rrms and oo 
yield good fits to the data, since fixing (3 = 1/Rrms gives x 2 /data point = 0.35, which is 
increased to only 1.02 for (3 = oo. 

2) In column 5 the use of only |.Fi| 2 as in ref. fU is clearly inferior with its x 2 /dp = 6.63. 

3) The normalisation factors needed to fit the experimental value of 0.63 /ib/sr at p v = 
0.246fm~ 1 are 0.87, 0.31, 1.50 for columns 3, 4, 5, respectively. This shows that I-F2I 2 is 1.7 
times stronger than |-F\| 2 and so could account for a significant part of the factor of 2.5 by 



which the model of Ref. [[1J] underestimated the experimental data. 

It should be added that there is a strong correlation between A and (3, e.g. for (3 = 
2/Rrms the minimum x 2 / dp is still 0.35 but with A=0.97. The dependence on the parameter 
A is also weak as it is with (3. Therefore, the main dependence may be expected to arise 
from the input values of the elementary rjN scattering amplitude. 

Unfortunately, the refinement in going from |-F\| 2 to I-F2I 2 gives less benefits with the 
other potential options. 

a) For potential III with (3 = 00 a x 2 /dp minimum of 0.58 occurs at A = 0.38 to be compared 
with x 2 /dp=0.71 for |-Fi| 2 i.e. little is gained by the refinement - in both cases a good fit 
being achieved to the data. Again there is a strong correlation between A and (3 with the 
above x 2 /dp=0.58 arising also for f3 = 1/Rrms and A = 0.14. 

b) Potential I gives already a good fit to the data using | -F\ | 2 with x 2 /dp=0.61. This cannot 
be matched by I-F2I 2 , which gives x 2 /dp=16 with (3 = 1/Rrms and A = 1. This only improves 
as (3 increases and A decreases i.e. finally back to \Fi\ 2 . 

c) Potential II is the worst combination. Here |-Fi| 2 gives x 2 /dp=6.5. In comparison I-F2I 2 
using (3 = 1/Rrms and A = 1 gives x 2 /dp=49, i.e. neither model gives a reasonable fit to 
the data. As with potential I, this only improves as the \Fi\ 2 limit is approached. 

The corresponding results with potential IV for 4 He are shown in Fig. 2. There it is seen 
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that |-Fi| 2 from Eq. fll8|) gives a visually better fit to the data and that ji^A = 0.88,(3 = 
^-/Rrms)\ 2 appears to produce too much energy dependence. However, it should be noted 
that here the experimental data have large error bars and exist only at a few energies. In 
the opinion of the authors, this should not be considered a fatal problem. Clearly some 
reduction of the experimental errors would be welcome to make these data more selective. 

So far the values of a(rjN) used are those suggested by experiment. However, these 
differ considerably amongst themselves with a(r/iV) = [0.3-0.6]+i[0.3-0.4]fm being a more 
reasonable estimate (see caption of Table |). In view of this, it is of interest to make a 
global variation of the input a(j]N) to recognize the optimal regions to fit the pd — > r/ 3 He 
cross section data within different model regimes. Such a calculation was performed for 
(3 = l/i?RMS an d A = 1, i.e. in a crude model where these are not varied. Fig. 3 shows 



the results for yx 2 /dp in the complex a(r]N) plane. In the hatched regions this parameter 
is smaller than unity and other contours show the values 2,3,. ..,10. It can be seen that 
there is a systematic change due to each correction introduced in this work into the optical 
model results - with all these additional effects being in the same direction. There may 
be a common area around a(r)N) = 0.4 + i0.3 fm for the optical model || without a Born 
background introduced in Eq. flT9D and for the full model, but elsewhere the models are 
exclusive. The 3 He data would allow in each model a valley of minimum x 2 i n different 
regions for a(rjN). So it is clear that (even assuming that the production mechanism were 
known) these data cannot uniquely determine the scattering length, although they set a 
strong constraint. It may be noted that similar fits could be attempted for the 4 He data. 
However, there the quoted experimental errors are so large that as such the fit would be 
useless. Even so, the energy independence of the production amplitude indicated by the 
four existing data points close to the dd — > 7/ 4 He threshold in very suggestive. It was 
not possible to produce this feature with any reasonable value of the elementary scattering 
lengths allowed by the above considered models for 3 He. Similar energy dependences in the 
4 He case are also obtained by Wilkin in Ref. [[01 . 
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IV. CONCLUSIONS 



This paper is in two distinct parts. In the first, the basic ?7-nucleon scattering length 
a(r]N) is converted into effective 77— 3,4 He scattering lengths a(r] 3 ' 4 He), which, in the second 
part, are then used to calculate the final state interactions in the pd — > 7/ 3 He and dd — > 7/ 4 He 
reactions. 

The step from a(rjN) to a(r/ 3 ' 4 He) is made in two stages using a multiple scattering 
expansion, the accuracy of which was first checked in the optical model limit - a limit that 
could be calculated directly from the Schrodinger equation (see Table |). Both the first 
stage, in which the replacement A — > (A — 1) is made, and the second stage, in which the 
scattering from a nucleon that is bound is taken into account, give large corrections that 
tend to go in the same direction. The overall effect is to give a(^ 3 ' 4 He)'s that are very 
different from those expected using the pure optical model (see Tables |H] and [II]). However, 
it should be added that this calculation ignores the effect of the possible presence of a sizable 
effective range in the basic r)N interaction. 

When the above a(^ 3,4 He)'s are used to extract the effect of final state interactions from 
the pd — > ?7 3 He reaction, it is found that only one (option IV) of the a(r/A/")'s proposed in 
the caption of Table | is able to give a good fit to the 3 He data - but not the less restrictive 
4 He data (see Figs. 1 and 2). 

In an attempt to see if there exist other values of a{r]N) that can give a good fit to 
3 He and, in addition, give a better fit to the 4 He data, a search was made in the region 
1.0 > Rea(r)N) > —1.0 and 1.0 > Im a(r]N) > 0.0. However, this did not produce any 
a(rjN) significantly better than the earlier option IV. If one may disregard the 4 He data 
either as too inaccurate or arising from too complex a reaction, it seems that the 3 He results 
indicate some potential for constraining the elementary rjN scattering length. 
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APPENDIX A: ISOSPIN STATES 

The 77-deuteron and ?7- 4 He systems are special cases. These are isospin systems. There- 
fore, some decay modes to the pion-nucleon channels are not allowed by isospin conservation. 
In the multiple scattering expansion this blocking is due to a cancellation of pionic waves 
emitted from several coherent sources. This effect has been shown to be important in co- 
herent ^-production processes [|17j . 



Here we first summarise briefly a two channel description of 77-nucleon scattering. Then 
we discuss the question of blocking the pionic channel in isospin systems. We follow the 
standard description in terms of a separable matrix T or V dominated by coupling to 
the iV*(1535) resonance. Let V be 



v " - A- < A1 > 



where M is the bare mass of the iV* and the fi are couplings to the different channels. The 
latter are functions of the channel momenta q v and q n . The scattering matrix T follows from 
the Lippman-Schwinger equation 

f = V + V(E - H + ie)- l f. (A2) 



The separability of the interaction (|A1|) then allows for the simple solution . 

T-— (AD 

il3 E-M -E k <G k >> [A6) 
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where 



<G k > 



(A4) 



E-E k ( q y 



The value of the total Re < G > yields an energy shift for the N*, while Im < G > 
determines its width. With a relativistically invariant normalisation N(q) for both the r\ 
and iV in Eq. (|A4|) one obtains for the partial width into channel k 



and the T m (0) = —a(j]N). The parameters of the coupling strengths and form factor 
ranges may be fitted to //-photoproduction (electro-production) and 7r-nucleon scattering 
data as well as to the N* decay properties. Analyses of this sort have been performed by 
several groups [|H,|9|,|T6ll ■ These differ slightly in the treatment of relativistic effects and on 
the (uncertain) input, and there is significant variation in the actual predictions for the 
scattering lengths a(rjN). 

Now, let us consider rj scattering on a correlated S = 1, T = pair of nucleons forming 
a quasideuteron state. The intermediate states in the pionic channels have T = 1 and 
so, due to the Pauli principle, the intermediate nucleons must be antisymmetric in space 
coordinates. This may reduce the available phase space and so lead to a blocking of virtual 
(or real) r] — it transitions. As a consequence the effective Im a(rjN) may be reduced in a 
nuclear medium. To allow for this effect we calculate the correction to the < G v > of Eq. 
(|A4|) due to this Pauli effect. An average quantity 



is calculated with an antisymmetrised free NN propagator - the average being taken over 
the AW ground state. Further, in this estimate a zero-range interaction is assumed between 
the meson and nucleons, which are considered to be fixed. In this way a correction term 
< AG n > is obtained in the form 



-j- = -Im < G k > 



q k {E)M N fl{q k {E)) 



(A5) 



AnE 



7 < U'i + / 2 )GW(/i + h) > 



(A6) 




(A7) 
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where 



A(q) = J dvi(t> 2 NN {n) 



2sin 2 (\^)-l 



-P(q). (A8) 



Here 4>nn is the initial NN wavefunction and p is the Fourier transform of the related 
density. For large systems this correction disappears, since < 2 sin 2 ( 3 ^ i ) >^ 1. But it 
could be sizable, if the inverse -Rrms of the system is comparable to the momenta involved. 
However, for Im < AG W >oc q-Kpilw) /^(ov) with q w w 2fm~ 1 one finds only a few per cent 
change of the N* width in the deuteron and in helium. This is so small a correction - also 
obtained at high momentum, where the wave functions tend to be uncertain - that it is 
reasonable to neglect its effect. 
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FIGURES 

FIG. 1. The pd — ► ry 3 He amplitude square |/(expt)| 2 defined in Eq. (21) plotted against the 



7/ momentum in the cm. system for the four elementary rjN amplitudes I-IV given in Table Q. 
Dashed curve: optical model; Dotted: optical model corrected by A — > A — 1 but described by Eq. 
(|i~8|); Dash-dot: off-shell effect also included in Eq. (|l~S|); Solid: the full model with corrections to 
the optical model and with the background term in the wave function (19). The data are from 
Refs. §|,|13|. 

FIG. 2. The dd — > ry 4 He amplitude squared for the elementary t]N amplitude IV. Dots: data 
from H (normalized as the total cross section). Curves as in Fig. 1. 

FIG. 3. Contour diagrams of \/x' 2 /dp for different models: a) optical model, b) optical model 
corrected by A — > A — 1 but described by Eq. ( [lq) , c) off-shell effect also included, d) the full 
model with corrections to the optical model and with the background term in the wave function 
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TABLES 

TABLE I. Comparison for rj-He scattering lengths o(r^He) (in fm) from various stages of the 
series expansion in Eq. (^) with the results from direct calculation with the corresponding optical 
potential [jl^]. The numbers are for 3 He and those in the brackets refer to 4 He. The results are 
illustrated with four sets of the rjN input: a(r]N) = 0.476 + i 0.279 fm (I), 0.579 + i 0.399 fm (II), 
0.430 + i0.394fm (III) and 0.291 + i 0.360fm (IV). In all cases R RM s = 1.788 [1.618] fm. 



a(r]N) 


a(IA) 


Q = R = 


R = 


a(r]A) 


a(opt) 


I 


1.89+il.ll 


-2.00+i3.01 


-1.94+i2.65 


-1.89+i2.60 


-1.87+i2.59 




[2.63+U.54] 


[-2.46+H.27] 


[-2.14+H.12] 


[-2.05+11.13] 


[-2.01+H.16] 


II 


2.30+U.59 


-2.41+U.94 


-2.16+H.71 


-2.08+H.70 


-2.06+il.72 




[3.20+i2.20] 


[-2.24+i0.83] 


[-1.91+10.79] 


[-1.81+10.83] 


[-1.79+i0.90] 


III 


1.71+U.57 


-1.67+12.12 


-1.56+H.93 


-1.52+H.92 


-1.51+H.93 




[2.38+i2.18] 


[-2.03+H.13] 


[-1.78+H.06] 


[-1.70+11.08] 


[-1.70+U.12] 


IV 


1.16+H.43 


-0.93+H.92 


-0.90+U.81 


-0.89+U.80 


-0.88+H.80 




[1.61+11.99] 


[-1.62+11.38] 


[-1.47+H.28] 


[-1.42+il.30] 


[-1.42+H.31] 
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TABLE II. The notation is the same as Table | except that the A n factors in P, Q and R are 
now replaced by (A — l) n giving aA-\{r)A) as indicated by Eq. (g). 



a(rjN) 


Qi = R' = 


R' = 


a A— i (t?j4) 


a(opt) 


I 


53+i4 27 


30+i4 19 


28+14 16 


-1 87+12 59 




[-3.01+i2.94] 


[-2.76+i2.55] 


[-2.67+i2.51] 


[-2.01+U.16] 


II 


-1.52+i4.40 


-1.61+i4.06 


-1.59+14.01 


-2.06+H.72 




[-3.03+H.84] 


[-2.67+il.64] 


[-2.56+H.65] 


[-1.79+i0.90] 


III 


-0.53+i3.35 


-0.61+i3.21 


-0.60+i3.18 


-1.51+11.93 




[-2.38+i2.23] 


[-2.16+i2.02] 


[-2.09+i2.02] 


[-1.70+11.12] 


IV 


-0.13+i2.36 


-0.16+i2.31 


-0.16+i2.30 


-O.88+1I.8O 




[-1.53+i2.25] 


[-1.45+i2.10] 


[-1.41+i2.09] 


[-1.42+11.31] 



TABLE III. The notation is the same as in Table [0| except that a(rjN) is now calculated 
off-shell using E sep + E mc = 19 [33] MeV corresponding to iK v = 0.581 [0.766] fin -1 . 



a(r]N, 0) 


a(r]N, off) 


R' = 


aA-i{r]A,oS) 


aA-i(r)A,Q) 


I 


0.39+i0.17 


2.01+i2.85 


1.99+i2.86 


0.28+i4.16 




[0.37+MU5] 


[-1.56+15.30] 


[-1.59+15.19] 


[-2.67+12.51] 


II 


0.47+i0.22 


1.36+i4.38 


1.32+i4.37 


-1.59+14.01 




[0.44+i0.18] 


[-3.14+13.88] 


[-3.07+i3.77] 


[-2.56+11.65] 


III 


0.39+i0.24 


0.93+i3.08 


0.92+i3.07 


-0.60+i3.18 




[0.37+i0.21] 


[-1.76+13.69] 


[-1.73+i3.62] 


[-2.09+i2.02] 


IV 


0.29+i0.26 


0.59+i2.17 


0.58+i2.17 


-0.16+i2.30 




[0.29+i0.23] 


[-0.78+i2.95] 


[-0.78+i2.93] 


[-1.41+12.09] 
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TABLE IV. The final state interaction factors |Fj(A,/3)| 2 in units of fib/sr for the elementary 
amplitude IV (afaJV) = 0.291 + i0.394), with A = 0.88 and /3 = l/i?RMS or oo. 



_Pr, l/(expt)| 2 \Ftf \F 2 ({3 = oo)\ 2 \Ftf 

0.051 0.53(0.02) 1.37 1.32 1.21 

0.115 1.07(0.03) 1.06 1.02 0.95 

0.166 0.86(0.015) 0.86 0.84 0.80 

0.202 0.74(0.014) 0.75 0.74 0.72 

0.246 0.63(0.020) 0.63 0.63 0.63 

0.295 0.50(0.016) 0.52 0.53 0.55 

0.337 0.45(0.018) 0.44 0.46 0.49 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9502022vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9502022vl 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9502022vl 



